Thermodynamic Properties of the Spin-1/2 Antiferromagnetic ladder 
Cu2(C2Hi2N2)2Cl4 under Magnetic Field 
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Specific heat (Cv) measurements in the spin-1/2 Cu2(C2Hi2N2)2Cl4 system under a magnetic field 
up to = 8.25r are reported and compared to the results of numerical calculations based on 
the 2-leg antiferromagnetic Heisenberg ladder. While the temperature dependences of both the 
susceptibility and the low field specific heat are accurately reproduced by this model, deviations 
are observed below the critical field Hc\ at which the spin gap closes. In this Quantum High Field 
phase, the contribution of the low-energy quantum fluctuations are stronger than in the Heisenberg 
ladder model. We argue that this enhancement can be attributed to dynamical lattice fluctuations. 
Finally, we show that such a Heisenberg ladder, for H > Hci, is unstable, when coupled to the 3D 
lattice, against a lattice distortion. These results provide an alternative explanation for the observed 
low temperature {Tc ~ 0.5/S' - 0.8/S') phase (previously interpreted as a 3D magnetic ordering) as 
a new type of incommensurate gapped state. 

PACS: 75.10 Jm, 75.40.Mg, 75.50.Ee, 64.70.Kb 



Wide interest is currently devoted to "gapped" spin 
systems, both experimentally and theoretically. In one 
dimension, 5' = 1 Haldane and alternating spin chains 
provide good examples of such systems. When magnetic 
frustration is present, similar situations can be found not 
only in one-dimension (ID) - the Ji — J2 model for in- 
stance - but also in two (2D) or three dimensions. An 
intermediate situation between ID and 2D is provided 
by the so-called "ladder" systems, which couple an even 
number of quantum spin {S — 1/2) chains. As for alter- 
nating and frustrated spin chains, the energy diagram is 
characterized by an energy gap A5 between the S* = 
ground state (GS) and the first magnetic S* = 1 excited 
state leading to characteristic magnetic and thermody- 
namic properties at low enough temperature, T < A5. 

The application of a magnetic field yields drastic 
changes in the energy spectrum. In particular, as a result 
of the Zeeman splitting undergone by the S = \ excited 
state, a second-order transition occurs at the critical field 
Hci — As/giJ-B, where g is the gyromagnetic ratio and 
fiB the Bohr magneton. Above Hqi, the GS becomes 
magnetic and a continuum of excitations develops giv- 
ing rise to "incommensurate" zero-energy fluctuations. 

Experimentally, the study of such a Quantum High 
Field (QHF) phase, i.e. for fl" > Hci, requires to work 
on systems having a relatively smallgap A5. Indeed, in- 
sulating ladders such as SrCu203 whose structure is 
closely related to the parent 2-dimensional cuprate anti- 
ferromagnets typically have spin gaps larger than lOOK. 
This explains that such a phase has rarely been investi- 



gated. As shown by recent studies, an interesting oppor- 
tunity is provided by the compound Cu2(C2Hi2N2)2Cl4 
(also known as CuHpCl) which is thought to behave as 
an ideal 2-leg spin-1/2 ladder system with a critical field 
of Hci — 7. sr. Magnetic measurements have been used 
first to characterize the magnetic parameters of the spin 
system. The behavior of the magnetic susceptibility, re- 
produced from Refs. in Fig. 0, is consistent with 
a gap of the order of A5 ~ IIK. Specific heat (Cy) 
measurements in CuHpCl have recently been performed 
under a magnetic field of up to 9T 1^. In low field 
{H < Hci), a single maximum is observed at relatively 
high temperature (T ~ J±), and, due to the presence of 
the energy gap, Cy decreases exponentially at low tem- 
perature [~ exp(— A5/T)]. In addition, a second order 
transition was shown to occur at very low temperature 
{0.5K <Tc < 0.8K) and was interpreted as the onset of 
3-dimensional (3D) magnetic order. 

In the present work, new specific heat measurements in 
a field (up to 8.25r) are presented which mainly focus on 
the QHF phase. They were performed by using a scan- 
ning adiabatic method. A small known power is applied 
to a high purity silicon sample holder and the temper- 
ature difference between the sample and a surrounding 
radiation screen is measured by a gold-iron thermocou- 
ple using a DC squid device as a current amplifier. A 
feedback network maintains the radiation screen at the 
same temperature as the sample, then strongly reducing 
the heat exchange process. The temperature rises slowly 
from 0.1 K to 8K, at a speed as slow as 10 mK/min. The 
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measurement of tlie temperature of the radiation screen 
then allows the specific heat of the sample to be calcu- 
lated. Such a slow drift rate in temperature ensures that 
all parts of the sample are in thermal equilibrium, and 
unlike pulsed methods the specific heat of non metallic 
materials with poor thermal diffusivities can be accu- 
rately measured. In the present work, four single crys- 
tals glued onto a mica plate (as in previous susceptibility 
experiments [3,4]) were measured. Each crystal weighed 
approximately 0.5 mg. The contributions of all the ad- 
denda - the mica , varnish, vacuum grease, and the silicon 
sample holder - were estimated and subtracted. The be- 
haviors observed are directly compared to the results of 
a numerical investigation based on the Heisenberg ladder 
model. In the QHF phase, a second maximum develops 
at low temperature. Above Tc, we observed deviations 
from the isolated ladder model (the contribution of the 
low energy " incommensurate" spin fluctuations occurs at 
larger fields) which, we argue, can be due to dynamical 
lattice fluctuations. We suggest an alternative explana- 
tion for the low temperature phase in terms of a new 
incommensurate gapped state. A calculation based on 
a Heisenberg ladder coupled to 3D (classical) phonons 
strongly supports this scenario. 
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FIG. 1. Theoretical fit of the temperature dependence of 
the magnetic susceptibility. The experimental data are taken 
from Refs. The anisotropy ratios Va = J±/J\\ as well as 

the magnetic coupling along the chain J — are indicated on 
the plots, (a) and (b) correspond to various set of parameters 
producing two different spin gaps As. 

The Hamiltonian we shall use to describe the com- 
pound is the Heisenberg model on a ladder deflned by 



(1) 



/3j 



where /3 (=1,2) labels the two legs of the ladder (oriented 
along the a:-axis), j is a rung index {j—l,...,L) and J|| and 
J_L are the bond strengths along and between the chains 



respectively. An applied field H in the Z-direction leads 
to an additional Zeeman term, Hz — -^gi^sH j 
with an average value g ~ 2.08 Q]. 

Our numerical approach is based on Exact Diagonal- 
isation techniques. At T = 0, clusters with size up to 
2 X 14 can be handled with the Lanczos algorithm allow- 
ing, after a proper finite size scaling procedure, for accu- 
rate determinations of the various physical quantities . 
At finite temperature, a full diagonalization of 2 x 6, 2 x 8 
and 2 X 10 ladders has been performed. According to 
previous literature, the anisotropy ratio Tq — J±/J|| lies 
around 5. In this regime, the spin correlation length is 
smaller than the system sizes so that finite size correc- 
tions become negligable. 

In order to test the choice of parameters in (|^), let 
us briefly consider the temperature dependence of the 
magnetic susceptibility Q. A comparison between the 
numerical and the experimental data is shown in Fig. ^. 
In fact, the quality of the fit is not very sensitive to the 
anisotropy ratio (in the range 4.75 < ra < 5.5). Pa- 
rameters producing a gap of Ag ~ 10. 8^ (as reported 
in Ref. Q) or As ~ 11.3 (as reported in Ref. |^) give 
excellent fits of the experimental data. As shown numer- 
ically or in a strong coupling approximation Q the 
experimental behavior of the magnetization vs H is well 
reproduced by a similar set of parameters. 

Concerning the specific heat measurements shown in 
Figs. ||, the exponential behavior at low temperature 
characteristic of the spin gap is suppressed at moderate 
magnetic fields. Above 7.5T, a broad maximum in Cy (T) 
builds up signaling the emergence of new low energy fluc- 
tuations [ro|. The numerical calculations of Cv(T) in 
Fig. H(b) based on the above ladder model (with param- 
eters leading to Hci — 7.7T) reveal qualitatively the 
same behavior. At low fleld, up to H = 6T, the agree- 
ment with experiment is very good, hence establishing 
the relevance of the ladder model (^ in this regime. In 
the QHF, however, the maximum observed in the theo- 
retical calculation appears at higher magnetic flelds than 
in experiment. Deviations from the theoretical behavior 
appear at low temperature for fields above 7.5T after the 
closing of the ladder gap. We argue here that this effect 
can be due to lattice fluctuations. In fact, it has been 
shown Jll| for uniform Heisenberg chains, in the context 
of spin-Peierls transitions, that an underlying spin-lattice 
coupling can lead to signiflcant deviations e.g. in the 
magnetic susceptibility which can be accounted for by 
an effective exchange coupling. As shown in Fig. ||(c), a 
behavior qualitatively similar to the experimental obser- 
vations can be obtained by using renormalized exchange 



couplings J'lf 



cJfi (a* =-L, II), c < 1, which, according 



to Ref. ||11| , is consistent with the effect of a coupling to 
the lattice. For increasing fleld above 7T, the renormal- 
ization parameter c decreases signaling an increasing role 
played by the lattice coupling. 

Motivated by the above discussion, we reasonably as- 
sume the presence of a magneto-elastic coupling along 
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FIG. 2. Specific lieat vs T for difFerent values of the mag- 
netic field, (a) Experimental data. (b,c) Comparison be- 
tween experiment (symbols) and theory (lines) using ra ~ 5, 
Jii = 2.63K. 



the chain direction by replacing the second term of 
Eq. (0) by 

^11 - ^11 + ^/3,,)S„/3 • S,+i,0 + ^kY^SI, , (2) 



/3,i 



where the second term corresponds to the (3D) lattice 
elastic energy and where the set of parameters {Sfj.j} 
(proportional to the atomic displacements) have to be 
determined by minimizing the total energy. 

Hamiltonian (||) can lead to a lattice distortion (i.e. 
that {Sp,j} =/= 0) in the strong coupling limit i.e. when 
Jj_ ^ J||. In this case, for H > Hci, by retaining 
the = 1 and — states only on each rung 
(see e.g. |l^,|l3|), the spin ladder reduces to a ID spinless 
fermion model with a hopping amplitude t = •^||/2 
and a nearest neighbor repulsion V — >/||/2. Physically, 
a particle corresponds, in the original spin language, to a 
Sz = 1 rung triplet excitation so that the effective band 
filling is directly proportional to the relative magnetiza- 
tion M/Msat in such a way that 2kF = 27r(M/Msat). 
If one neglects the short range repulsion V between the 
particles (this should be justified for low particle den- 
sity i.e. small magnetization), as in the usual Peierls 
transition, a modulation of the hopping amplitude t = 
J||/2 of wavevector 2kF and magnitude S opens a gap 
at the chemical potential and leads to an energy gain 
AE cx (5^ In (const /(5), for S ^ 1. For arbitrary large K, 
the minimum of the total energy is then obtained for an 
equilibrium value S ~ exp (— constif). 

In order to explore the relevance of the previous sce- 
nario, a numerical investigation on finite clusters is re- 
quired. The following study has been restricted to T = 
and simple ratios for M/Msat hke 1/3, 1/2 or 2/3. The 
minimization of the total energy using expression (^ can 
be realized by an Exact Diagonalization technique sup- 
plemented by a self-consistent procedure ||l^. A typical 
GS configuration is shown in Fig. ^(a) for a magneti- 
zation M = I A/sat- The lowest energy is obtained for 
a perfectly symmetric modulation of the two chains i.e. 
<5/3j = Sj. The variation of the distortion Sj along the 
ladder is correlated with that of the spin density and is 
a periodic function of period Xp — Mgat/M. For the 
special case M/M^at = 1/2, the distortion becomes com- 
mensurate and corresponds to a simple dimerization of 
the lattice, similar to the D-phase observed in spin-Peierls 
chains such as CuGeOs in the absence of a magnetic field. 

In order to study the stability of these modulated 
phases we have defined a critical elastic constant as 
Kc = lini\Sj\^Q{AE / J2j <^j}i where AE corresponds to 
the magnetic energy gain due to the equilibrium dis- 
tortion pattern. The distorted phase is then stable for 
K < Kc- Our results displayed in Fig. |^(b) as a func- 
tion of the ratio J^^l J\\ show that Kc increases with sys- 
tem size; although it is difficult to extrapolate our results 
to the thermodynamic limit, they clearly establish that 
a small spin-lattice coupling leads to modulated struc- 
tures. It is interesting to notice that, for a small M/Msat 



3 




•—•2x8, M=1/2 
■—■2x12, M=1/2 
0-0 2x6, M=1/3 
^^2x12, M=1/3 
+ - + 2x6, M=2/3 
*- + 2x12, M=2/3 



10 11 12 



(b) 



— ^ — * - ^ + 




3 4 

r. 



FIG. 3. (a) Equilibrium modulation 1+Sj and correspond- 
ing average spin density {{Sf^i + Sf_2)) vs the position along 
the ladder direction calculated on a 2 x 12 cluster for Va ~ 5, 
K = 0.6 and M/Msat = 2/3. (b) Critical value Kc as a func- 
tion of the anisotropy ratio ra — J±/J\\ for M/Msat ~ 1/3, 
1/2 and 2/3. The system sizes are indicated on the plot. 



(i.e. for H just above Hci), a simple 2kp modulation 
is expected (in this regime, the short range repulsion V 
becomes irrelevant) while, for M/Msat close to 1/2 more 
complicated incommensurate structures of "soliton lat- 
tice" type (i.e. involving an infinite number of higher 
harmonics) |16| should be stabilized. All these incom- 
mensurate structures bear strong similarities with the I- 
phase of the spin-Peierls systems and have similar 
properties as e.g. the existence of a gap fl^ . 

To conclude, we have shown that the specific heat data 
for H < Hci are quantitatively well described by an 
isolated Heisenberg ladder model. Deviations from the 
predictions of this model observed for H > Hci and 
low temperatures are attributed to the effect of a small 
magneto-elastic coupling to the 3D lattice. We have 
shown numerically that, in this regime, the Heisenberg 
ladder becomes unstable against a lattice distortion lead- 
ing to a new gapped incommensurate phase. 
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